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Abstract 

We calculate the unintegrated gluon distribution at low intrinsic transverse momenta 
and its parameters are found from the best description of the SPS and LHC data on 
the pp collision in the soft kinematical region. It allows us to study the saturation 
of the gluon density at low Q 2 more carefully and find the saturation scale. 



1. Introduction 

As is well known, hard processes involving incoming protons, such as deep-inelastic 
lepton-proton scattering (DIS), are described using the scale-dependent parton den- 
sity functions. Usually, these quantities are calculated as a function of the Bjorken 
variable x and the square of the four-momentum transfer q 2 = —Q 2 within the 
framework of popular collinear QCD factorization based on the DGLAP evolution 
equations [1]. However, for semi- inclusive processes (such as inclusive jet production 
in DIS, electroweak boson production [2, 3], etc.) at high energies it is more appro- 
priate to use the parton distributions unintegrated over the transverse momentum 
k t in the framework of fc r factorization QCD approach [4], see, for example, reviews 
[5, 6] for more information. The k t - factorization formalism is based on the BFKL 
[7] or CCFM [8] evolution equations and provides solid theoretical grounds for the 
effects of initial gluon radiation and intrinsic parton transverse momentum k t . The 
theoretical analysis of the unintegrated quark q(x, k t ) distribution (u.q.d.) and gluon 
g(x, k t ) distribution (u.g.d.) can be found, for example, in [9]- [13]. In this paper we 
estimate the u.g.d. at low intrinsic transverse momenta kt < 1.5 — 1.6 GeV/c and its 
parameters extracted from the best description of the LHC data at low transverse 
momenta pt of the produced hadrons. We also show that our u.g.d. similar to the 
u.g.d. obtained in [9, 10] at large k t and different from it at low k t . The u.g.d. is 
directly related to the dipole-nucleon cross section within the model proposed in [9] , 
see also [13]- [18], that is saturated at low Q or large transverse distances r ~ 1/Q 
between quark q and antiquark q in the qq dipole created from the splitting of the 
virtual photon 7* in the ep DIS. So, we find also a new parameterization for this 
dipole-nucleon cross section, as a function of r, from the modified u.g.d. and analyze 
the saturation effect for the gluon density. 
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2. Inclusive spectra of hadrons in pp collisions 



2..1 Unintegrated gluon distributions 

As was mentioned above, the unintegrated gluon density in a proton are a subject of 
intensive studies, and various approaches to investigate these quantities have been 
proposed [7]. At asymptotically large energies (or very small x) the theoretically 
correct description is given by the BFKL evolution equation [7] where the leading 
ln(l/x) contributions are taken into account in all orders. Another approach, valid 
for both small and large x, is given by the CCFM gluon evolution equation [8]. 
It introduces angular ordering of emissions to correctly treat the gluon coherence 
effects. In the limit of asymptotic high energies, it almost equivalent to BFKL [7], 
but also similar to the DGLAP evolution for large x ~ 1. The resulting unintegrated 
gluon distribution depends on two scales, the additional scale q is a variable related 
to the maximum angle allowed in the emission and plays the role of the evolution 
scale ji in the collinear parton densities. Also it is possible to obtain the two-scale 
involved unintegrated quark and gluon densities from the conventional ones using the 
Kimber-Martin-Ryskin (KMR) prescription [5, 6]. In this way the kr dependence 
in the unintegrated parton distributions enters only in last step of the evolution, 
and usual DGLAP evolution equations can be used up to this step. Such procedure 
is expected to include the main part of the collinear higher-order QCD corrections. 
Finally, a simple parameterization of the unintegrated gluon density was obtained 
within the color-dipole approach in [9] on the assumption of a saturation of the 
gluon density at low Q 2 which successfully described both inclusive and diffracting 
EPA scattering. This gluon density xg(x, k 2 ,Ql) is given by [9, 10] 

xg(x, k t , Q ) = 3a ° R 2 k 2 exp (-R 2 (x)k 2 ) ; R = ^- (—) , (1) 

where a = 29.12 mb, a s = 0.2, Q = 1 GeV, A = 0.277 and x = 4.1 • 10~ 5 . This 
simple expression corresponds to the Gaussian form for the effective dipole cross 
section a(x, r) as a function of x and the relative transverse separation r of the 
qq pair [9]. In fact, this form could be more complicated. In this paper we study 
this point and try to find the parameterization for xg(x, k t , Qo), which is related to 
a(x, r), from the best description of the inclusive spectra of charge hadrons produced 
in pp collisions at LHC energies and mid-rapidity region. 

2. .2 Quark-gluon string model (QGSM) including gluons 

As is well known, the soft hadron production in pp collisions at not large trans- 
fer can be analyzed within the soft QCD models, namely, the quark-gluon string 
model (QGSM) [19]-[22] or the dual parton model (DPM) [23]. The cut n-pomeron 
graphs calculated within these models result in a reasonable contribution at small 
but nonzero rapidities. However, it has been shown recently [24, 25] that there are 
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some difficulties in using the QGSM to analyze the inclusive spectra in pp collisions 
in the mid-rapidity region and at the initial energies above the ISR one. However, 
it is due to the according to Abramovsky-Gribov-Kancheli cutting rules (AGK) [26] 
at mid-rapidity (y ~ 0), when only one-pomeron Mueller-Kancheli diagrams con- 
tribute to the inclusive spectrum ph(y ~ 0,p t ). To overcome these difficulties it 
was assumed in [24] that there are soft gluons or the so called intrinsic gluons in 
the proton [27], which are split into qq pairs and should vanish at the zero intrinsic 
transverse momentum (k t ~ 0). The total spectrum p h (y ~ 0,p t ) was split into 
two parts, the quark contribution p q (y ~ 0,p t ) and the gluon one and their energy 
dependence was calculated [24, 25] 

p(Pt) = p q (x = 0,p t ) + p g (x = 0,p t ) = g(s/s ) A (f) q (0,p t ) + (g(s/s A - a nd ) (f> g (0,p t ) . (2) 
Here 

p q (x = 0, p t ) = g(s/ s ) A (f) q (0, p t ) (3) 

and 

p g (x = 0,p t ) = (g(s/s A - ow) 4>g(0,p t ) (4) 

The following parameterization for <j) q (0,p t ) and <j)g(0,pt) was found [24]: 

g (O,Pt) = A q exp(-b q p t ) 
<f>g(Q,Pt) = Ag^p~ t exp(-b g pt), (5) 

where so = lGeV 2 ,g = 21mb,A = 0.12. The parameters are fixed from the fit to 
the data on the p t distribution of charged particles at y = [24]: A q = 4.78 ± 0.16 
(GeV/c)" 2 , b q = 7.24 ± 0.11 (Gev/cT 1 and A g = 1.42 ± 0.05 (GeV/c)" 2 ; b g = 
3.46 ± 0.02 (GeV/c) _1 . Figure 1 illustrates the best fit of the inclusive spectrum of 
charged hadrons produced in pp collisions at ^Js =7 TeV and the central rapidity 
region at the hadron transverse momenta p t < 1.6 Gev/c; the solid line corresponds 
to the quark contribution p q , the dashed line is the gluon contribution p g , and the 
dotted curve is the sum of these contributions ph given by Eq.(2). 

In Figs. (2,3) the inclusive spectra of the charged hadrons produced in pp colli- 
sions at the mid-rapidity region and yfs =2.36 GeV, 7 Tev, 500 GeV, 900 GeV are 
presented with the inclusion of the calculations within the perturbative LO QCD, 
see details in [25]. One can see from Figs. (1-3) that the calculations within the soft 
QCD including both the quark contribution (Eq.(3)) and the gluon one (Eq.(4)) 
and the contribution, which corresponds to the LO QCD calculation, results in the 
satisfactorily description of these spectra in the wight region of the initial energies. 
It stimulates us to study the form of the gluon contribution 4> g (x 0,ptt) related to 
p g (x 0,pt) in detail and find the information on the distribution of the soft gluons 
at small transverse momenta. 
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Figure 1: The inclusive spectrum of the charged hadrons as a function of pt (GeV/c) 
in the central rapidity region (y = 0) at =7 TeV at pt < 1.6 GeV/c compared 
with the CMS [28] which are very close to the ATLAS data [29]. The solid line is 
the quark contribution p q (x = 0,p t ) (Eq.(3), the long dashed curve corresponds to 
the gluon one p g {x = 0,p t ) (Eq.(4), the dotted line is the sum of the quark and 
gluon contributions (Eq.(2). 

2.. 3 Modified unintegrated gluon distributions 

As can be seen in Figs. (1-3) the contribution to the inclusive spectrum at y ~ due 
to the intrinsic gluons is sizable at low p t < 2 GeV/c, e.g., in the soft kinematical re- 
gion. Therefore, we can estimate this contribution within the nonperturbative QCD 
model, similar to the QGSM [19]. We calculate the gluon contribution 4> g (x ~ 0,p t ) 
entering into Eq.(5) as the cut graph (Fig. 4, right) of the one-pomeron exchange 
in the gluon-gluon interaction (Fig. 4, left) using the splitting of the gluons into the 
qq pair. The right diagram of Fig. 4 corresponds to the creation of two colorless 
strings between the quark/ ant iquark (q/q) and antiquark/quark (q/q). Then, after 
their brake qq are produced and fragmented to the hadron h. Actually, the calcula- 
tion can be made in a way similar to the calculation of the sea quark contribution 
to the inclusive spectrum within the QGSM [19], e.g., the contribution (p g (0,p t ) is 
presented as the sum of the product of two convolution functions 

4>g(x,p t ) = F q (x + ,p ht )Fq(x-,p ht ) + Fq(x + ,p ht )F q (x-,p ht ) , (6) 
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Figure 2: The inclusive spectrum of charged hadron as a function of pt (GeV/c) 
in the central rapidity region (y = 0) at yfs =2.36 TeV (left) and y/s =7 TeV 
(right) compared with the CMS [28] and ATLAS [29] data. The solid line (soft 
QCD (quarks)) is the quark contribution p q (x = 0,p t ) (Eq.(3)), the dotted curve 
(soft QCD (gluons)) corresponds to the gluon one p g (x = 0,p t ) (Eq.(4), the long 
dashed line (SQCD(quarks+gluons)) is the sum of the quark and gluon contributions 
(Eq.(2), the short dashed curve (Perturbative QCD) corresponds to the perturbative 
LO QCD [25] and the dash-dotted line (SQCD+PQCD) is the sum of the calculations 
within the soft QCD including the gluon contribution (Eq.(2)) and the perturbative 
LO QCD. 



where the function F q ( q )(x +) pht) corresponds to the production of final hadrons from 
decay of qq string. It is calculated as the following convolution: 

Fq(q)(x±,Pt,Pht) = J ' ± dxi J d 2 k lt f q{q) (xi,k lt )G q{ q)^ h (j^>Pto ~ k t)^j , ( 7 ) 

Here G q ^ h (z,k t ) = zD q{q) ^ h (z,k t ), D q(q) ^ h (z,k t ) is the fragmentation function 
(FF) of the quark (antiquark) to a hadron h, z — x±/xi,k t = pnt — kt, x± = 

0.5(^Jx 2 + x\ ±x),x t = 2^J (ml + pi)/ s. The distribution of sea quarks (antiquark) 
f q (q) is related to the splitting function V g -> qq of gluons to qq by 

fq(q)(z,h) = I g(zi,k t ,Q )Vg^ q q( — )—, (8) 
Jz Z\ Z\ 

where g{z\, k u , Qq) is the u.g.d.. The gluon splitting function V g ^ qq was calculated 
within the Born approximation. In Eq.(8) we assumed the collinear splitting of the 
intrinsic gluon to the qq pair because values of k t are not zero but small. 

Calculating the diagram of Fig. (4) (right) by the use of Eqs.(4-8) for the gluon 
contribution p g we took the FF to charged hadrons, pions, kaons, and pp pairs 
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Figure 3: The inclusive spectrum of charged hadron as a function of pt (GeV/c) 
in the central rapidity region (y = 0) at yfs =540 GeV (left) and y/s =900 GeV 
(right) compared with the UAl [30] and ATLAS [29] data. The solid line (soft QCD 
(quarks)) is the quark contribution p q {x = 0,p t ) (Eq.(3), the dotted curve (soft QCD 
(gluons)) corresponds to the gluon one p g [x = 0,p t ) (Eq.(4), the long dashed line 
(SQCD(quarks+gluons)) is the sum of the quark and gluon contributions (Eq.(2), 
the short dashed curve (Perturbative QCD) corresponds to the perturbative LO 
QCD [25] and the dash-dotted line (SQCD+PQCD) is the sum of the calculations 
within the soft QCD including the gluon contribution (Eq.(2)) and the perturbative 
LO QCD. 
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Figure 4: The one-pomeron exchange graph between two gluons in the elastic pp 
scattering (left) and the cut one-pomeron due to the creation of two colorless strings 
between quarks /antiquarks (right) [19]. 



obtained within the QGSM [31]. From the best description of p g [x ~ 0,Pht, see hs 
parameterization given by Eq.(5), we found the form for the xg(x, k t , Qo) which was 
fitted in the following form: 

xg(x, kt, Qo) = . 2 3ff ° . Ci(l - xfv x 
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R 2 (x)kf + C 2 (R (x)k t ) a ) exp {-R (x)k t - d(R (x)k 
The coefficient G\ was found from the following normalization: 

rQl 



9(x,Ql) 



dk 2 g(x,k 2 ,Ql) , 



(9) 



(10) 



and the parameters 



0.7; G 



2 - 2.3; A = 0.22; b g 



12;d = 0.2;C 3 = 0.3295 



were found from the best fit of the LHC and SPS data on the inclusive spectrum of 
charged hadrons produced in pp collisions and in the mid-rapidity region at p t <1.6 
GeV/c, see the dashed lines (SQCD (quark+gluons)) in Figs. (1-3) and Eq.(5). 
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Figure 5: The unintegrated gluon distribution xg(x,k t ,Q )/Co as a function of k t 
at x = xo and Q = l.GeV/c. The dashed curve corresponds to the original GBW 
[9, 10], Eq.(l), and the solid line is the modified u.g.d. given by Eq.(9). 

Figure 5 presents the modified u.g.d. obtained by calculating the cut one- 
pomeron graph of Fig. 4 and the original GBW u.g.d. [9, 10] as a function of 
the transverse gluon momentum k t . Here Co = 3a / (47r 2 a s (Qo)) . One can see that 
the modified u.g.d. (the solid line in Fig. 5) is different from the original GBW 
u.g.d. [9, 10] at kt < 1.5 GeV/c and coincides with it at larger k t . This is due to 
the sizable contribution of p g (Eqs.(4,5)) to the inclusive spectrum p(pt) of charged 
hadrons produced in pp collisions at the LHC and SPS energies in the mid-rapidity 
region, see the dashed lines (soft QCD(gluons)) in Figs. (1-3). 

Let us also note that, as is shown recently in [32], the modified GBW given 
by Eq.(9) does not contradict the HERA data on the longitudinal structure func- 
tion Fi(Q 2 ) at low x, the charm structure function F2 C (x,Q 2 ) and the bottom one 
F 2b (x,Q 2 ). 
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3. Saturation dynamics 

According to [9], the u.g.d. can be related to the cross section a(x, r) of the qq dipole 
with the nucleon. This dipole is created from the split of the virtual exchanged 
photon 7* to qq pair in the e — p deep inelastic scattering (DIS).This relation at the 
fixed Ql is the following: [9]: 

a(x,r) = J d ^{l-Urk t )}xg{x,k t ) (11) 

Inserting the simple form for xg(x,k t ) given by Eq.(l) to Eq.(ll) on can get the 
following form for the dipole cross section: 

* G B W (x,r) = * {l-exp(- T ^} (12) 

However, the modified u.g.d. given by Eq.(9), inserted to Eq.(ll) results in the more 
complicated form for a(x, r): 

r) = {l - ex P (-^ - ^ } , (13) 

where b x = 0.045, b 2 = 0.3. 

There are a few forms for the dipole cross sections suggested in [13]-[18]. The 
dipole cross section can be presented in the general form: [9]: 

a(x,r) = cr g(f 2 ) , (14) 

where f = r/(2Ro(x)). The function g(r 2 ) was presented in the form [13] 

g (f 2 ) = f 2 ln (l + (15) 

or in the form [17]: 

g (r 2 ) = 1-exp j-f 2 Zn(^ + e^j (16) 

which both of them are saturated when r grows. The function g(f 2 ) was also pre- 
sented in the form of type [14, 15]: 

g(f 2 ) = ln(l +f 2 ) (17) 

that is not saturated when r increases. Figure (6) illustrates the dipole cross sections 
ct/oo at x = xq which are saturated at r > 0.6 fm, see [13, 16, 17]. They are 
compared with our calculations (solid line, Modified a) given by Eq.(13). The solid 
line in Fig. 6 corresponds to the modified u.g.d. given by Eq.(9) the application 
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Figure 6: The dipole cross section a/a at x = Xq as a function of r. The green 
line "GBW cr" is calculation of [9], the red line "Modified a" corresponds to our 
calculation; the curve "Nikolaev,Zakharov a" is the calculation of [13]; the line 
"Albacete,Marquet a" is the calculation of [17]. 



of that allowed us to describe satisfactorily the LHC data on inclusive spectra of 
hadrons produced in the mid-rapidity region of pp collision at low p t . Therefore, the 
form of the dipole-nucleon cross sections presented in Fig. 6 can be verified by the 
last LHC data on hadron spectra in soft kinematical region. 

Comparing the solid line (Modified a) and dashed curve (GBW a) in Fig. 6 one 
can see that &Modif(x,r) given by Eq.(13) saturates faster than aGBw(%, r ) given by 
Eq.(12) when the qq dipole distance r increases. If Ro = l/GeV(x/xo) x ^ 2 , according 
to [9], then for saturation scale has the form Q s ~ l/-Ro — Qso{ x o/ X ) 7 where 
Q s q = 1 GeV=0.2 fin -1 . The saturation in the dipole cross section (Eq.(12) sets in 
when r ~ 2R or Q s ~ (Q s0 /2)(x /x) x ^ 2 ). Comparing the saturation properties of 
the Modified a and GBW a presented in Fig. 6 one can get slightly larger value for 
Q s q in comparison to Q s0 = 1 GeV/c. 



4. Conclusion 

Assuming creation of soft intrinsic gluons in the proton at low transverse momenta k t 
and calculating the cut one-pomeron graph between two gluons in colliding protons 
we found the form for the unintegrated gluon distribution (modified u.g.d) as a 
function of x and k t at the fixed value of Qq. The parameters of this u.g.d. were 
found from the best description of the SPS and LHC data on the inclusive spectra 
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of the charged hadrons produced in the mid-rapidity pp collisions at low p t . It was 
shown that the modified u.g.d. is different from the original GBW u.g.d. obtained 
in [9, 10] at k t < 1.6 GeV/c and it coincides with the GBW u.g.d. at k t > 1.6GeV/c. 

Using the modified u.g.d. we calculated the qq dipole-nucleon cross section aModif 
as a function of the transverse distance r between q and q in the dipole and found 
that it saturates faster than the oqbw obtained within the GBW dipole model [9]. 
It allowed us to find the saturation scale Q s for the gluon density that is larger than 
the one obtained in [9]. Moreover, we showed that the satisfactory description of 
the LHC data using the modified u.g.d. can verify the form of the dipole-nucleon 
cross section and the property of the saturation of the gluon density at low Q 2 . 
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